Simulating 2+ Id Lattice QED with dynamical matter using ultracold atoms 
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We suggest a method to simulate lattice compact Quantum Electrodynamics (cQED) using ultracold 
atoms in optical lattices, which includes dynamical Dirac fermions in 2+1 dimensions. This allows 
to test dynamical effects of confinement as well as 2d flux loops deformations and breaking, and to 
observe Wilson-loop area-law. 



Recent progress in the manipulation and control of 
atomic systems has boosted the interest in the emerg- 
ing field of quantum simulation [T]-[3]. So far, most of 
the theoretical and experimental effort has concentrated 
in simulating condensed matter systems. The field of 
quantum simulation, however, may also have a strong 
impact in high energy physics (HEP) [2], lending us the 
possibility of observing intriguing phenomena emerging 
from the standard model, or even allowing us to answer 
some questions which cannot be addressed with standard 
techniques of lattice QCD. In contrast to its condensed 
matter physics counterpart, the field of quantum simula- 
tion of HEP models is almost unexplored, and presents 
new theoretical and experimental challenges. 

In the standard model of HEP, forces are mediated 
through gauge bosons, which are the excitations of gauge 
fields. Thus, gauge fields, either abelian or non-abelian, 
play a central role in particle physics. As such, gauge the- 
ories have been a subject of continuous research over the 
last decades, theoretically, numerically and experimen- 
tally. One important property of such theories is known 
as confinement [SHE], which is manifested in real- world 
QCD, in the structure of hadrons: Free quarks cannot be 
found in nature, due to the phenomenon of confinement, 
which "holds them together" forming hadrons. Confine- 
ment is known to hold in non-abelian gauge theories (not 
only QCD), but it is also manifested in abelian theories, 
such as U(l) - QED: in the lattice case - Compact QED 
(cQED). There, in 3+ld, a phase transition takes place 
between a confining phase (for strong couplings) and a 
Coulomb phase (for weak ones), and in 2+ld confine- 
ment takes place for all values of the coupling constant 
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Recently, some suggestions for simulations of QFT and 
HEP models have been proposed. This includes the simu- 
lations of dynamic scalar [12] (vaccum entanglement) and 
fermionic fields [13] (the interacting Thirring and Gross- 
Neveu models). Simulation proposals for fermions in Lat- 
tice QFT, either free or in external non-dynamical gauge 
fields include Axions and Wilson fermions [14], Dirac 
fermions in curved spacetime |15j and general quantum 
simulators of QFTs and topological insulators [TB]. On 
the other hand, some simulations were proposed for dy- 



namic gauge fields but with no fermions, using BECs |17j 
or single atoms [TS] in optical lattices, where the first 
simulates the abelian Kogut-Susskind Hamiltonian [10] 
and the latter - a truncated ("Spin-Gauge") version of it 
(other examples of truncated models are p~9l42 1| ) . More 
simulations of pure-gauge U(l) theories with ultracold 
atoms have been recently suggested [251 123] • 

The next step is the inclusion of dynamical matter 
(fermions) in the model, allowing for simulation of full 
cQED. This is of special interest, both for the possibil- 
ity of probing confinement in a dynamic matter theory, 
as well as for the exploration of problems which are not 
amenable of a numerical description, due to the well- 
known Grassmanian integration sign problem [24J. The 
first proposals for the simulation of a fermionic system 
with a dynamic gauge field [Ml [25] are, so far, restricted 
to 1+ld, where the emerging physical phenomena are 
limited due to the absence of magnetic fields and of mul- 
tiple paths connecting two points (see, however, the QED 
(non-lattice) work [26 ). 

In this paper, we propose a spin-gauge model contain- 
ing dynamic fermions, that allows for simulation of non- 
trivial gauge field dynamics for 2+1 cQED, including 2d 
spatial effects such as flux loops. We show how to con- 
struct a gauge invariant model of a U(l) gauge field cou- 
pled to 'naive' fermions, discuss how to create various 
initial states, and suggest several possible experiments, 
such as the observation of broken flux lines and manifes- 
tation of Wilson's area- law of confinement [5J, realizing 
the gedanken experiment proposed in |27j . 

The Simulating System. Consider a 2d spatial optical 
lattice [3j, filled with fermionic and bosonic atoms. Each 
vertex n = (m, 71%) is occupied by a single fermion, either 
of species C or D, described by the local Hilbert spaces 
spanned by eigenstates of the fermionic number operators 
N^' D (with the annihilation operators c n , d n ). The two 
species form a spinor at each vertex: 

We define the local charge [28J, 

On = V4V^n - 1 = N° + ~ 1 (2) 
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Each link (in directions k = 1,2) is occupied by a single 
boson, of one of 21 + 1 species, forming an angular mo- 
mentum multiplet, described by eigenstates of L\ n (See 
also the figure in [2"5]V 

We call this system "the primitive theory" since the 
gauge- invariant effective theory will arise from it, by 
adding a constraint, as will be shown later, but first we 
shall describe the (separate) primitive theories for the 
bosons and the fermions. Similar fermionic theories can 
be found, for example, in [HlU6| . 

The bosons primitive dynamics is governed by the 
Hamiltonian 



h: 



E 

n.k 



(id) 



where the second sum is on nearest-neighbor (intersect- 
ing) links. 

The fermions primitive dynamics is governed by the 
lattice Dirac Hamiltonian 

H = ivJ2 (V>W n+ £ - H.C.) + Mj>t ^ ( 4 ) 

n.fe n 

Imposing Gauge invariance. We wish to impose gauge 
invariance to the system. In order to do so, we con- 
strain the generators of gauge symmetry G n to zero, by 
adding the "Gauss's Hamiltonian" Hq to the two primi- 
tive Hamiltonians. This Hamiltonian involves an interac- 
tion between the two primitive systems. The generators 
of local gauge transformations are 



G n — L\ +L 2 Z +L 1 



-L 



,-(-1 



ni+n 2 



Qn (5) 



and the constraint is implemented by the Gauss Hamil- 
tonian, Hq = A^G^ and thus A must be the highest 

n 

energy scale: X ^> (i, 0, fi, rj, M. 

As A is the largest energy scale, we wish to derive 
an effective theory for Hq's ground sector [30] . which 
will introduce gauge invariance to the system. In the 
first order, we get the bosonic "electric Hamiltonian" 
He = v-'zn) > an d the fermionic "Mass Hamilto- 

nian" H M = Mj^ipt^ipn 
j i 

In the second order we get four contributions: (1) 
The bosonic "magnetic" plaquette terms: Hb = 
_20i£ ( L i +tuL *_ Ll 2 iZ n + H.C.) as well as the 

bosonic H' B , which has no equivalent in cQED, but is 
gauge-invariant and contributes no interesting dynamics 
(see [TH] for details). (2) The "minimal coupling" Dirac 
terms: H D = H^ n+{c L k s>n - H.C.) (where 

n,k 

s = (— )™ 1+ ™ 2 ). The other two contributions are less im- 
portant, and are described in the supplementary material 



Next we make the sign change L 



k 

z,n 



Besides inverting the sign of L z on odd links, we also 



swap the L± operators there (leaving L x invariant). That 
results, as in the pure spin-gauge theory, with correct 
signs in the plaquette term 

H B = (^n^n+i^n+An + B.C) (6) 

n 

and the gauge generators 

Gn = (— 1)" 1+ " 2 \ Lg tn + L 2 z n — L l z n _j — l? z n 2 - Qn 

(7) 

This also gives rise to the correct (minimally coupled) 
Dirac Hamiltonian 



(3) 



H n = 



it]j3 



(8) 



where we identify the Dirac matrices as the Pauli 
matrices a\ 2 (The Dirac j3 matrix is a z ). 

Thus, we have obtained the effective Hamiltonian (ne- 
glecting H' B ): 



H 
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where a is a rescaling factor: a = -^fi = 

M 



_ il 2 (l + l) 2 n 2 g 2 _ 
~ A — 

277^-y/i(i+i ) _ ^ introducing the gauge theory coupling 
constant g and fermions mass m [53]. Equation ([9| is 
the Kogut-Susskind Hamiltonian of lattice cQED with 
dynamical Dirac Fermions. 

We wish to emphasize that the previous effective 
Hamiltonian calculation does not depend on I, and thus 
the primitive theory and the constraint are achieved for 
other values of I, the same effective Hamiltonian results. 

Realization of the model. In order to realize the model 
for 1 = 1, for the intra-bosonic terms one can use the 
methods presented in [IB], which generalize [3T]. Almost 
all the required intra-bosonic terms in ^ (including the 
ones in Hq, but not the linear L x terms) are obtained if 
one sets, in [18] terms, q n = everywhere. To get the 
linear L x terms of ([3]) , one must add an on-site species ro- 
tation term which can be achieved using Raman lasers, 

E y/l(l + l)-m(m + l)al 



(where 



=-/ 



a m is the local mth bosonic species annihilation opera- 
tor). This does not affect the effective calculation since 
these are local terms. 

For details on the fermionic interactions, see [29]. The 
main issue is now the realization of the Boson-Fermion 
interaction terms in Hq. The required terms, for each 
vertex n, are of the form 



L, 



2.11—1 2 . n — 2 , 

(10) 

where £ (n) = — 2(— 1)" 1+ " 2 For a start, we choose the 
atomic hyperfine levels of the fermionic degrees of free- 
dom: C and D fermions belong to the hyperfine lev- 
els \Fq — ^,mc) , I-Fd = §, Wo) respectively, where 



3 



mc,niD are constant on each vertex: The hyperfine lev- 
els are arranged such that the other values of m, on each 
hyperfine manifolds, are of a much larger energy, and thus 
are not accessible if they are not initially populated. 



n, m c = {-l) ni+ " 2 \m c \ ,m D = (-l) r " ' \m D \, we 
are only left with one type of term, which is, due to the 
third condition, -2A (-1)™ 1+ " 2 (n£ + N£) L z , which is 
the requested interaction. 

No Interactions vacuum and excited states. Sup- 
pose first that the bosonic and fermionic interactions 
(Hb,Hd) are turned off - i.e., the system is pure-gauge 
(the charges are static and uncoupled) and in the strong 
coupling limit. Hence the gauge field vacuum does not 
contain any excited links: all the bosonic atoms are in 
their m — state. This is the exact QED vacuum in 
the extreme strong limit, since no plaquette terms con- 
tribute. 

As for the fermionic content, each vertex can con- 
tain four different fermionic contents. The state with 
no charges corresponds to the "Dirac Sea" - i.e., all the 
vertices are filled with D atoms. In such a state, the 
local "mass" is —M. Since energies are measured above 
the Dirac sea, we would like to measure these "masses" 
above —M, i.e. shift the energies with the product of M 
with the total number of vertices. This is the vacuum in 
case of no charges. The three excited states correspond 
to different simulated fermionic contents, as summarized 
in the following table; Our "particle" is q, with a positive 
charge of +1, and its "anti-particle" is q, with a negative 
charge of —1. Note that the masses are relative to — M: 



\n 1 +n 2 



The ( 10 ) terms are derived as scattering interactions 
between two particles: bosons with Fj, = 1 and each 
of the fermionic species. The first quantized interaction 
Hamiltonian corresponding to such scattering processes 
is H sc = J29fPf> where the summation is on the total 

F 

angular momentum of the two scattered particles, {<?f} 
depend on the S-wave scattering length, and thus are 
tunable using optical Feshbach Resonances |32ti34] . and 
Pf are projection operators to the subspaces of different 
total angular momenta. These operators can be built 
using the different values of Fc ■ F b , Fb ■ Ft, for each 
value of total F, and thus one can express the scattering 
Hamiltonians for C, D as 



H sc — Co + CiFc ■ F b 

H? B = D Q + D{F D ■ F b 



D 2 (F D -F b y (11) 



where d = d (gi ,53 J , D, = Di (gi,ga,gs 

These Hamiltonians are next plugged into the sec- 
ond quantized interaction terms / d 3 x^ i '^Hf c ^ 
where $z are the fermionic and bosonic second quan- 
tized wavefunctions. Due to the energy spectrum of 
the fermions, the only possible processes are such with 
i = k = mc- Thus no angular momentum transfer can 
take place for the bosons as well, and hence j = I = m b , 
and only the z components of the angular momentum 
vectors contribute. We finally get, on each vertex, 

m c m b ) cl 1 c a a'l lb a mb + 
m D m b + D 2 f {m D ,m b )) d^d 

(12) 

where the bosonic operators correspond to each of the 
neighboring links, Ci,Di are products of the coeffi- 
cients Ci,Di and the appropriate overlap integrals, and 
/ (m£>,m b ) is quadratic. 

In order to get (10 1, we impose conditions on the 



Co [91,91 

c( 



\m D \ 



scattering coefficients gF (controlled by Feshbach res- 
onances): D 2 (gi,ga,g §) = 0, Dq [ g ± . a ±. <r_L 
and L>i (51,33,55 

\mc\ = — 2 A the first one eliminates 
the quadratic terms. The second one results in terms of 
the form C (Ng + Nj?) N b . However, the bosonic ef- 
fective Hamiltonian sets that N b = 1 everywhere around 
the bosonic lattice, and thus this term is merely the total 
number of fermions in the system, which is an ignorable 
constant in the Hamiltonian. Setting that at the vertex 
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Initial state preparation. The system is initially pre- 
pared in the gauge + fermions ground state, which is an 
eigenstate of Hq+He+Hm (Bosons - m = everywhere, 
and all the vertices are occupied only with D fermions). 
Then, these three Hamiltonians can be turned on, with- 
out changing the state but imposing the constraint. 

There are several interesting initial states one can ob- 
tain, even before turning Hf> on. These include: (a) Ob- 
taining the QED static vacuum (up to first order, if one 
works with I = 1), by increasing Q adiabatically, keeping 
20 j (i+i) ^ ^. Obtaining a "loop sea" by increasing 
O adiabatically again, but with 20 1 ^> ^1, dressing 

to a state with a lot of loops, but with larger amplitudes, 
and then, by lowering the value of J7, one can obtain 
strong regime QED dynamics again; (c) Creating, using 
single-addressing lasers [351 13S]; zeroth order excitations: 
one can create large closed flux-loops as in the pure gauge 
case, but also and/or place C charges and empty vertices 
to form zeroth-order flux tubes. Afterwards Hb can be 
turned on adiabatically to obtain QED dynamics. 

Inclusion of Matter Dynamics. Given any of the initial 
states (a),(b) or (c), the field-matter interaction Hb can 
be turned on (in order to remain in the strong limit, 
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FIG. 1. Examples of gauge field and matter dynamics, as 
explained in the text, for an initial state (using single ad- 
dressing lasers) of a single large flux loop (a), (b-f) show 
possible outcomes of first order processes: (b-d) result from 
single operations of Hb ■ (b) loop decreasing, (c) loop increas- 
ing, (d) exciting a non-connected plaquette; (e-f) result from 
single operations of Hd '■ (e) flux loop breaking and creation 
of two charges, (f) creation of a meson outside the loop with 
a new flux tube connecting them. 



the parameters must obey 



A A 

|29|). Then, an interplay between the fermion masses and 
the fermionic interaction terms introduces changes in the 
flux loops/tubes structure, if one waits long enough. For 
example, suppose 2M = fiL, for some integer L. Then, 
if there is a flux tube/loop longer than L, we expect it 
to break to two flux tubes, generating two new fermions, 
transferring the energy of the " disappearing links" to the 
masses of the new fermions. On the other hand, if the 
fermion mass is small (or zero) we expect the flux tube 
to break and form new fermions constantly. 

Another possibility is to move slowly (externally) one 
of the charges of an initially prepared flux-tube and 
stretching it, which should lead to it breaking apart. 

Examples of matter dynamic processes are shown in 
figure 0. 

Wilson's Area law observation. Another possible mea- 
surement which can be done in this system is a test of 
Wilson- loop area law, as a probe for confinement. This 
can be done by an interference of two mesons in superpo- 
sition. This cannot be done in real-world QED, where it 
is only a gedanken experiment |27j . Here we propose to 
realize it in the quantum simulation. For that, we need 
the fermions to be static and their mass is unnecessary, 
and hence HM,Hf are turned off (and thus, effectively, 
Ho as well). In fact, for a proof of principle of this ex- 
periment, Hb is not needed too. 

We start, in the extreme strong limit, with a state of a 
single flux-tube with length R over the Dirac Sea - \R), 
emanating between the vertices (m,n) and (m + R, n), 
such that Q mtn = l,Q m +R t n = — 1. This state is just a 
meson in the strong coupling limit, and it is an eigenstate 
of the Hamiltonian. Then, as in [27] . we wish create 
a superposition of two states, \R) and \R+ 1): i.e., to 
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FIG. 2. Description of the area-law experiment in space-time, 
as described in the text. The interference phase depends on 
the green area. 



transfer the C fermion, with a probability ~, from (to, n) 
to (to — 1, n). We shall treat this subspace as a two level 
system, denoted by \\) = \R),\\) = with the 

Hamiltonian H = /i \f) 

The required superposition is generated as in Ramsey 
intcrfcrometry [29J : we apply on the initial state \l) the 
unitary operation U — e~ l ^ ay , and get the state U \l) = 
— |4-x)- Next, wait a long time r>T, during which the 
state evolves to ^e~ itiRT (||) - e^ T Operating 

on the state with U again yields (|J, X ) + er l ^ T |f x )j 
(neglecting the global phase). Switching back to the a z 
basis, recalling the definitions of our effective two-level 
system, we get that the probability to find a fermion C 
in the vertices (to, n) , (to — 1, n) are 



2 ( uLT\ 
P m ,n = sin I — — I 



i Mn-l,n 



= COS 



(13) 



since A = LT (above L = 1), this manifests the area 
law in a confining phase, and thus, within this system 
a realization of the gedanken experiment proposed there 
is possible. In future generalizations, this method may 
serve as phase transition probe. 

Summary. In this paper, we have proposed a method 
to add dynamic fermions to the Spin-Gauge model intro- 
duced in [TB]. This allows to simulate 2+ld cQED using 
ultracold atoms in optical lattices, enabling the obser- 
vation of confinement in a dynamic matter theory, and 
suggesting a detour to the well known Grassmanian inte- 
gration sign problem, encountered in Monte-Carlo simu- 
lations [23] . Simulations of our model on a small lattice 
may be used to check the effects of imperfections of the 
parameters on the dynamics. We have suggested several 
possible measurements, which can show the emergence of 
dynamic charges, the breaking of flux-tubes and the area 
law behavior of the confining potential. 
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SUPPLEMENTAL MATERIAL 
Tailoring the fermion interactions 




FIG. 3. The simulating lattice. Every link is occupied by a 
single boson (yellow circle) and every vertex by either zero, 
one or two fermions (blue circles). 

The intra-fermionic terms (in the fermionic primitive 
Hamiltonian) can be obtained if one tailors the fermionic 
Raman transitions and chemical potentials properly. 
Note that the required hopping terms, achieved using 
Raman lasers, include "species rotation", i.e.: for exam- 
ple, when a C fermion hops to a neighboring vertex it 
becomes a D - this is manifested by the Pauli matrices 
a Xl cr y in the interactions. However, this rotation can be 
avoided, if one initially defines slightly different primi- 
tive theories as we shall next describe. Other required 
intra-fermionic terms are the ones on each vertex, result- 
ing from AQ^ in Hq\ some add up to a constant - the 
total number of fermions in the system, and can be disre- 
garded, and the others are local fermion scattering terms: 
2\N C N D . 

v n n 

Now we describe an alternative method to get the 
same effective theory as in the paper, but with same- 
species fermionic tunneling and slightly different primi- 
tive Hamiltonians. We start by defining another spinor 
at each vertex, rather than ^> n , 



Xn = 




(14) 
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(note that the charges remains the same, as they are 
invariant for local fermionic species swapping.) and the 
primitive Hamiltonians are 
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n,fc n <i,j) 

where the second sum in on nearest-neighbor (intersecting) links, and 

H = (^nX„ + l - H.C.) + V J2 + B.C.) + (-1)" 1+ " 2 X ^Xn (16) 

n n ii 



if one now performs the transformation 

Xn = <7£ 1+ " 2 V>n 



(17) 



and repeats the effective calculation for this model, the 
same effective Hamiltonian will result. 

Since the charges remain the same, one could produce 
the bosonic-fermionic interactions exactly as described 
in the paper, where ip was the fundamental spinor, since 
they depend only on the invariant charge. 

More terms of the effective Hamiltonian 

In the paper we have described two types of second- 
order terms in the effective Hamiltonian. We hereby de- 
scribe other two second-order contributions, which are 
irrelevant for our purposes: 

1. Double operation of Hf (fcrmionic tunneling): 
Such terms may represent back-and-forth tunneling 
along a link. Some add up to the total number of 
fcrmions in the system, which is an ignorable con- 
stant of motion, while others arc C-D swap terms 

2 

along a link. These terms are of order ^- and thus 
can be attenuated (compared to the Dirac terms) 
if one sets 77 <C /?• 



2. Double operation of Hx- we get, on every 
link, a term of the form — -j^- L_} = 

— ^f- (L 2 — L 2 ). the first term is a constant. The 
second is just a renormalization of [i. 



Rescaling the Hamiltonian parameters 

In equation (9), we have defined the rescaling factor 
a which translates the experimental parameters into the 
simulated system parameters, introducing the coupling 
constant g and the fermion mass to. The relation between 
the two families of parameters is given, as written in the 
main text, by 



a 



2 _ 4l 2 {I + if tt 2 g 2 _ 2 V l3y / l (I + 1) _ M 

g 2 ^ 1 A A to 

(18) 

This relation contains the translations of the electric, 
magnetic, Dirac and mass terms. 



One should note that the electric, magnetic and mass 
part contain constants of the theory - g and to, and thus 
have some more freedom in their definitions. The Dirac 
kinetic part, however, contains, in the cQED Hamil- 
tonian (setting the lattice spacing a = 1) no physical 
constant, and thus the experimental parameters within 
it define a. Thus, if we wish to examine the strong 
coupling limit, for example, demanding g ^> 1 from 
the electric part, we get that » |f, which means 

11 » Wf"i 



netic part we get 40 1 



On the other hand, from the mag- 
•C a, which corresponds 

to 2^ + i)3 ^ wra This yiclds the strong limit 

condition, in the presence of dynamic fermions, which is 
stated in the paper: 



2n 2 l 2 (l + lf r)Py/l (T+Tj 



A 



A 



(19) 



One should also note that all this is true if one wishes 
to simulate "regular" cQED. However, if one chooses the 
parameters differently, the theory will still be gauge in- 
variant (as long as the scale hierarchy is kept and the 
effective expansion is valid) , but not the " regular" cQED 
theory. 



Comparison of simulating and simulated processes 

In terms of the simulating system, we shall denote the 
spinors as 



(20) 



Due to the Pauli matrices which couple nearest- 
neighbor vertices, it is only possible to either both create 
a particle and an anti-particle in the two ends of a link, or 
to both annihilate such a pair. This is phrased in terms 
of the simulating system as the fact that creation of a C 
fermion in a vertex means the annihilation of a neighbor- 
ing D, and vice-versa. Let us consider these equivalence 
of the simulating and simulated systems in some basic 
dynamic processes. 

Fig. 4 shows some of these processes - considering the 
initial and final states in both terms, and the parts of 
the Hamiltonian which are responsible for them (only 
the specific relevant parts of the spinors are included). 
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Creating the unitary transformation for the Ramsey 
interferometry 

Using Raman lasers, we introduce the interactions 

(^ c m,n c rn—l,n ~ c rn-l,n c m,nj (21) 

where rj <C A, which are similar to the ones in Hf. These 
interactions violate the constraint, and thus we can use 
them in order to construct effective Hamiltonian terms 

"^2\ (^ C rn,n c m-l,nL( m -l t n)- — cj n _ 1 „C m! „Z/( m _i jrl ) + ^ 

which can be rewritten as 



(Without loss of generality, we assume that (m — 1, n) is 
an even vertex) . We turn these terms on and off suddenly, 
keeping them open a period of time r = ^^ A . r 3> 
A -1 , which justifies the effective treatment. On the other 
hand, we choose the time T (see the main text) to satisfy 
t <C T, making H y strong and impulsive. Thus, the time 
evolution of the state when H y is on is equivalent to the 
unitary operation U = e~ l %' 7y , and if the initial state 
was ||), we now get U \i) = — |i x ). 



(22) 
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FIG. 4. Possible fermionic processes, both in terms of the 
simulated system - QED (the upper rows) and of the simu- 
lating system (the lower rows), including the relevant oper- 
ators. The initial states are on the left (a-c) and the final 
are on the right, (a) Creation of a particle-anti-particle flux 
tube (a meson), (b) Breaking a flux-tube, (c) "Stretching" 
a flux-tube. Notice that it is not really stretching, since a 
pair of particles appear in the middle, (d) Moving a me- 
son: first process (c) occurs, and then a second stage, shown 
in (d), takes place. In this figure, in terms of the simulat- 
ing system (lower rows), |0) c |0) D = |0), |0) o = \C), 
\Q) a \l) D = \D),\l) c \l) D = \CD). 



